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Accurate quartic anharmonic force fields for CF4 and SiF4 have been calculated using 
the CCSD(T) method and basis sets of spdf quality. Based on the ab initio force field 
with a minor empirical adjustment, the vibrational energy levels of these two molecules and 
their isotopomers are calculated by means of high order Canonical Van Vleck Perturbation 
Theory(CVPT) based on curvilinear coordinates. The calculated energies agree very well 
with the experimental data. The full quadratic force field of CF4 is further refined to the 
experimental data. The symmetrization of the Cartesian basis for any combination bands 
of Td group molecules is discussed using the circular promotion operator for the doubly 
degenerate modes, together with tabulated vector coupling coefficients. The extraction of 
the spectroscopic constants from our second order transformed Hamiltonian in curvilinear 
coordinates is discussed, and compared to a similar procedure in rectilinear coordinates. 

I. INTRODUCTION 

Carbon tetrafluoride (also know as Freon-14) and tetrafluorosilane are molecules of wide ranging industrial, 
environmental, and economical interest. While the use of carbon tetrafluoride is being phased out togethe* 
with that of other CFCs, it has a number of other important industrial and technological applications.!!] 
These include use as an isolator and extinguisher gas in high- voltage applications, and as an etching gas 
in the semiconductor industry.EJ The formation of CF4-. . . O2 adducts was considered as a model system for 
the use of liquid perfluorochemicals in artificial blood.13 Tetrafluorosilaae, SiF4, is a precursor (e.g. by glow 
discharge in SiF 4 -H2 mixtures) of amorphous Si-F-H semiconductors^ as well as for plasma deposition of 
low-dielectric Si-O-F thin solid films. El It is the main by-product of beam etching of both semiconductors with 
fluorine, and silicon dioxide with fluorocarbons.a In addition to being used to monitor the above processes 
[see e.g. Ref. ^|, SiF 4 has been proposedo as a remote volcano monitoring probe, since its presence can 
be measured using open-path FT-IR spectroscopy of f3(SiF 4 ). Finally, the isotopic separation of 30 Si (in 
30 SiF 4 ) by infrared multi-photon dissociation (IRMPD) of natural-abundance Si2F@ using CO2 lasers has 
been reported on a preparative scale.Q 

A key step for elucidating these processes is to understand the equilibrium structure, spectroscopy and 
energetics of these molecules. For these reasons the vibrational spectroscopy of the above molecules has 
been the subject of several theoretical and experimental studies. It is also the subject of the present paper 
in which we apply a combination of an accurate ab initio electronic structure treatment and a high-order 
CVPT treatment of the vibrational problem to the CF 4 and SiF 4 molecules. 

Early gas-phase spectroscopic measurements of the CF 4 moleculeBO have been reviewed and supplemented 
by Jones, Kennedy, and Ekberg (JKE).II9 Earlier, Jeannotte et a/.tll obtained spectra in liquid argon solution, 
which as expected differ considerably from the JKE ijesults due to solvent shifts. Early high-resolution work 
on the molecule was stimulated by the announcement£3 of mid-infrared CF 4 lasing based on the V2 + ^4 — > vi 
transition. Esherick et aZEj determined v\ to high resolution by inverse Raman spectroscopy. High-resolution 
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Raman spectroscopy of v-i was carried out by Lolckjll and of 2^2 a^fl v\ by Tabyaoui et aZ.il The v-i + v± 
combination band was measured to high-accuracy by Patterson et alS3 and by Poussigue,ei al.tJ The band 
was studied extensively by Jones et al.p by Tarrago et al.j^3 and by McDowell et aZ.,t2l the latter of whom 
determined v± for all three isotopic species with Doppler-limited resolution,, pp ppler-limited measurements 
of the 2v\ + combinatio.n-.band were published by two different— groups EBB A Fermi type _Linteraction 
exists between 2i/4 and ^3;E3 detailed analyses of the 1/3 and 2^30 and of the z/3/2^4 polyadEfl have been 
carried out. 

The spectroscopic studies of CF4 have been accompanied by several theoretical (saectroscopic) inves- 
tigations. Hawnonic valence fo^ee fields were first derived by by Duncan and MillsE3 and by Chalmers 
and McKean.Ej Jeannotte et alEU fitted a harppnic potential, supplemented by diagonal cubic and quartic 
stretching force constants only, to their earlieilij liquid argon solution measurements. BrodersenEa obtained 
a cubic force field from experimental rotation- vibration data. A quartic force field is said to be partly refined, 
to rovibrational levels below 140&-cm _1 , but detailed information about the fitting is unavailable so far.E^I 
At higher energies, Boujut et alEQ consider local and normal mode behavior in the z/3 and ^4 ladders of 
tetrahedral XY4 and octahedral XYg molecules, with CF4 being considered amons-the former. In their 
analysis they employ data from Ref. ^ as computerized into the STDS data bank.Eil The 712/3 ladder has 
also been found to be a good testing ground for comparing experiment and theory. Using their results of 
high-resolution proton energy loss spectroscopy, Maring_ef aZ.,E3 review earlier data fox-i^, 2^3, and 3^3. 
They also summarize earlier unpublished data by HeenanLj who derived two sets of Hechto-type anharmonic 
constants (to be denoted Heenan L-and Heenan II in the remainder of the paper) from fitting Urey-Bradley 
type force fields to the Jeannottetil and JKE data, respectively. Maring et al. propose some alterations 
based on their measurements, including an apparently exceptionally large value faz-X^ = —9.1 cm -1 . 

Early work on the spectroscopy of SiF4 has been reviewed by McDowell et aZ.,E3 who note that research 
on the molecule was stimulated by the fact that the 1/3 fundamental overlaps the P branch of the 9400 nm 
band of the CO2 laser. McDowell and coworkers reported a partial set of anharmonicittfpCpnstants as well. In 
recent theoretical work, the low lying vibrational spectrum of SiF4 has been modeledEjO using an algebraic 
approach. 

Patterson and PineE3 determined fBp=0. 13676(3) cm -1 , whence ro=l. 55982(17) AcB. In later work, Takami 
and KuzecJ and later Jorissen et a/.EJ substantially revised -Bo upward (to 0.137780439(92) cm" 1 , consistent 
with a substantially shorter r^—l. 55404 A). As this paper was being finalized, we received a preprint by 
Demaison et alEn, in which r e =l. 5524(8) is derived from a combination of ab initio and experimental results. 

In the present paper we re-investigate the spectroscopy of these molecules using a combination of accurate 
ab initio anharmonic force fields and advanced techniques-for solving the vibrational Schrodinger equation. 
In an earlier electronic structure study, Martin and TaylorcJ revised the heat of vaporization of silicon (which, 
among other things, is required for any ab initio or semiempirical calculation of the heat of formation of 
any silicon compound) from a benchmark ab initio calculation on SiF4 from a very precise fluorine, bomb 
calorimetric measurements of the heat of formation of SiF4(g) and a benchmark ab initio calculations of the 
total atomization energy of said molecule. In this paper, the focus is restricted to the molecular vibrations. 
In part, we were motivated to study-,these molecules due to the success of a recent study in which, an ab 
initio quartic force field for methaneEj obtained using large basis sets and coupled cluster methods served as 
the starting point for several spectral refinement studies using variational methods,!^ low-order perturbation 
theory,c3 and high-order canonical Van Vleck perturbation theory (CVPT).C3 In general, we find that only 
the quadratic force constants and perhaps the geometry needs to be refined, and most of the remainder of 
the force field can be constrained to the ab initio values. 



II. AB INITIO ANHARMONIC FORCE FIELD 

All electronic structure calculations were carried out using MOLPRO 970 running on an SGI Origin 2000 
minisupercomputer at the Weizmann Institute of Science. 

Electronic correlation was treated at the CCSD(T) level [coupled cluster with all single and double sub- 
stitutions CC,CSD)l£j supplemented with a quasiperturbative estimate o£-the contribution of connected triple 
excitationsEit] method, using the implementation thereof in MOLPRO.E3 This method is knowr£j to yield 
correlation energies that are very close to the exact basis set correlation energy (i.e. the full configuration 
interaction result in thB-same basis set) as long as correlation is predominantly dynamical in character. As 
pointed out previously^ correlation in SiF4 is essentially rpurely dynamical in character, and the same holds 
true for CF4. For the record, values for the T\ diagnosticEj (a measure for the importance of nondynamical 
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correlation) are only 0.0118 for SiF4 and 0.0119 for CF4, respectively. (In all calculations reported here, only 
valence electrons were correlated.) 

Geometry optimizations were carried out by univariate polynomial interpolation. At the reference geome- 
try, a quartic force field in curvilinear symmetry coordinates was determined by recursive application of the 
central finite difference formula with step size 0.01 A or radian. (Quartic contamination was removed from 
the quadratic force constants.) The symmetry coordinates are the same as those used in the previous study 
on CH4.EI 

The resulting curvilinear internal coordinate force field was then transformed to rectilinear normal coor- 
dinates, then a standard-, second-order rovibrational perturbation theoryEJ was carried out by means of the 
SPECTRO program .Pti 57 ! The alignment conventions for the anharmonic constants of a spherical top follow 
the-work of HechtLj and general formulae for these constants were taken from the paper by Hodgkinson et 
al.EIi (See Section III.B for a discussion about the definition of the v% normal coordinates.) The implementau. 
tion in SPECTRO was slightly modified to allow deperturbing for an arbitrary number of Fermi resonances,Ej 
and a simple routine implementing Hccht'sEj energy level formulas up to four vibrational quanta was added. 
Similar procedure has also been implemented in the curvilinear formalism, and the relationship between 
the two formalism will be discussed in Section III.B. As in past work (e.g. Ref. ^), correct alignment was 
verified by slightly (0.00001 a.m.u.) perturbing the masses of two atoms, then repeating the analysis in the 
asymmetric top formalism. 

Given the n 3 N 4 scaling (with n the number of electrons correlated and N the number of basis functions) of 
CPU time for a CCSD(T) calculation, the large number of valence electrons correlated (n—32) imposes certain 
restrictions on the size of the basis set for a complete quartic force field. .We-settled on the standard Dunning 
cc-pVTZ (correlation consistent polarized valence triple zeta) basis setElJ'EJ (150 contracted Gaussian- type 
functions, CGTFs) for CF 4 , while for SiF 4 , we opted for the cc-pVTZ+1 basis set (159 CGTFs), where 
the suffix "+1" denotes the addition of a single high-exponent d function on second-row elementsE3 to 
accommodate inner polarization effects. The latter has been repeatedly shown to be essential for accurate 
computed properties of second-row molecules with polar bonds£j l.3 and the addition of a mere five functions 
to the basis set increases CPU time insignificantly for our purposes. It should be pointed out that even with 
these relatively compact basis sets, the evaluation of each force field took the equivalent of six weeks of CPU 
time on the SGI Origin. 

For purposes of validation, we in addition calculated geometries and harmonic frequencies at the CCSD(T) 
level with a number of other basis sets. _Epr CF4, we considered the (diffuse-function) augmented cc-pVTZ 
basis set aug-cc-pVTZ of Kendall et alE3 (230 CGTFs), while for SiF4, we considered both the standard 
cc-pVTZ basis set (154 CGTFs) and the aug-cc-pVTZ+1 basis set (239 CGTFs), in which the "+1" suffix 
again denotes the addition of an inner-polarization d function. The results are collected in Tables | and ||. 

For CF4, the main deficiency appears to be that the asymmetric stretching frequency L03 ks-seriously over- 
estimated. This is a known problem for polar bonds, particularly those involving fluorine £3~E3 Increasing 
the basis-£|et size to cc-pVQZ (275 CGTFs) or even cc-pV5Z (455 CGTFs) as a rule does not remedy this 
problemE3'Ea — and given the ./V 4 scaling behavior would increase CPU time requirements by one and two or- 
ders of magnitude, respectively, and would therefore be essentially impossible. The use of anion basis set (i.e. 
aug-cc-pVnZ) on at least-tlie electronegative atoms themselves generally leads to a dramatic improvement in 
basis set convergence ,Ej~E3 and here too we see that CCSD(T)/aug-cc-pVTZ frequencies for CF4 are much 
closer to the final adjusted values than the CCSD(T)/cc-pVTZ values. .While in HF, this sensitivity extends 
to the anharmonicitiesj£3 in heavy-atom systems like FNO and ClNOEZl the anharmonicities are generally 
unaffected by this change. We have therefore chosen not to calculate the entire anharmonic force field at the 
CCSD(T)/aug-cc-pVTZ level, which would have increased CPU time requirements approximately fivefold. 

Interestingly, the effect of adding the anion functions is much less pronounced for SiF4 than it is for CF4. 
The effect of adding the inner polarization function is fairly pronounced on the bond distance (a decrease 
of almost 0.01 A), but effects on the harmonic frequencies are fairly minor. Here too, we judged that the 
CCSD(T)/cc-pVTZ+l level represented the best compromise between accuracy and computational cost for 
the entire anharmonic force field. 

Purely ab initio calculation of the geometry and the harmonic part of the potential would presumably be 
feasible as far as-.basis set convergence and inner-shell correlation corrections are concerned, but we know 
from experienced that imperfections in the CCSD(T) method might still account for (usually positive) errors 
of several cm -1 in the frequencies. Full configuration interaction calculations are out of the,question for this 
system, and hence we have followed a different track similar to previous work on silane,t3 and refined the 
'raw' CCSD(T)/cc-pVTZ and CCSD(T)/cc-pVTZ+l force fields in the following way. 

(1) Using standard second-order rovibrational perturbation theory, the four vibrational fundamentals for 
the most abundant isotopomer were obtained, together with the equilibrium and ground-state rotational 
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constants B e and B , respectively. Hence (as r e y/B e /B^) the computed r was obtained. 

(2) Using the internal coordinate force field given above, the four diagonal quadratic force constants and 
r e were then iteratively adjusted such that the fundamentals and ro obtained by second-order rovibrational 
perturbation theory coincide with the experimental values to the precision given. (.F47 and all cubic and 
quartic force constants were held at their ab initio co mpu ted values.) This process converged after three 
iterations, and the final force fields are given in Table III. This force field was then used as input for the 
higher-order vibrational analysis. 

As can be seen in Tables [I] and [n], the adjustment only entails relatively minor changes in the geometry 
and harmonic frequencies. In addition, overall agreement between the available experimental energy levels 
and the values computed using standard second-order perturbation theory is such that we conclude that 
the CCSD(T)/VTZ(+1) quartic force fields are of sufficient quality for our purposes. (Since the computed 
energy levels are essentially identical to those obtained by m eans of second-order CVPT if no resonances are 
accounted for, the reader is referred to Tables VI and VII] and the relevant sections below for details.) 

Concerning the equilibrium geometries, we iHpte that our final r e -walues for CF4 and SiF4 agree with 
the experimentally derived values of BrodersenEa and Demaison et al.^n respectively, to within their stated 
uncertainties. 



III. CALCULATION OF VIBRATION ENERGY LEVELS 



In this and the following two sections, in order to assess the accuracy of the ab initio force field, the 
vibrational energy levels of CF4 and SiF4 are calculated using high order canonical Van Vleck perturbation 
thcory(CVPT) and compared to the experimental data. The calculation is similar to that which has been 
recently applied to the methane and its isotopomers. Details about this calculation can be found in Ref. 
and only the major procedures will be summarized here. I— . 

Firstly the exact J = vibrational Hamiltonian including the pseudo-potential V' terrrJld is expanded 
in terms of the curvilinear normal coordinates based on the Morse coordinates for the stretch and angle 
extension coordinates for the bend. The redundancy relation between the six-bond angles is properly taken 
into account in the expansion. As in the case of methane and its isotopomers,c3 the kinetic energy operators 
are expanded to fourth order, even though CVPT is carried out to a higher order. We choose this scheme 
because the kinetic contributions of order higher than four are believed to be small, and the number of terms 
rapidly increases with the order of the expansion. The quartic potential which is already in expanded form 
is used throughout this work. It contributes to the CVPT up to the second order. It is re-expressed in the 
coordinates described above so that the new expansion of the potential agrees with the original potential up 
to quartic terms. 

The expanded Hamiltonian is separated in the form 

H v = + AffW + A 2 H (2) + ■ ■ • + \ n H^ n \ (1) 

where A is the perturbation parameter. A succession of canonical transformations of the form 

K v = exp{i\ n [S {n) , }} ■ ••exp{a 2 [S' (2 >, ]} exp{iA[5 (1) , }}H V (2) 

are then applied to the Hamiltonian of Eq_,|pt), where the are chosen such that the matrix representation 
of K v is block diagonal through order nOQ The blocks are characterized by one or more polyad quantum 
numbers 

N t = mjVj, (3) 

i 

where the integers mi define the block. These integers are chosen such that the resonance interactions are 
not perturbatively decoupled. Details regarding the choice of the polyad quantum numbers will be discussed 
later with respect to specific molecules. 

Two points, regarding our analysis of this Hamiltonian are now considered in more detail. First, in order to 
obtain and analyze the eigenvalues of the individual blocks we diagonalize each block in a symmetrized basis 
to reduce the matrix size and to help assign the calculated levels. Second, the transformed (effective) Hamil- 
tonian of Eq. (|^) at second order is particularly interesting as its coefficients are the quartic anharmonicity 
constants and resonance constants, which have been conventionally used by spectroscopists to highlight the 
features of the vibrational spectra. In the following two subsections we discuss the basis .svmmetrization and 
extraction of spectroscopic constants as they have not been given in our previous workJ!3 
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A. Symmetrization in the Cartesian basis for the Td point group 

There are generally two ways of choosing basis functions for the doubly and triply degenerated normal 
modes for the T4 group molecules. The first choice (angular momentum basis) are eigenfunctions of the 
vibrational angular momenta. The second choice (Cartesian basis) is the products of harmonic oscillator 
functions including different Cartesian components of the degenerate vibrations_In the angular momentum 
basis, the off-diagonal terms are the spherical tensor operators defined by HechtEj and their matrix elements 
are calculated using the Wigner-Eckart theorem. However it is more difficult to treat vibrational resonances 
in the angular momentum basis, and therefore the alternative Cartesian basis is more often used. Another 
advantage of the Cartesian basis is that the matrix elements arc readily calculated since they all reduce to 
one-dimensional harmonic oscillator matrix elements. For both basis, the symmetrization is nonobvious and 
we are going to discuss the symmetrization in the Cartesian basis below. 

The Cartesian basis denoted by \n-\_\ n2a, n2b] Ti3 X , n^y, n^ z \ n^ x , n^y, n^) will be symmetrized into the irre- 
ducible representations of the Td point group: Au-,~Aa, E, F\ and F2. Various equivalent approaches can be 
used, such as the projectiorLjoperator technique JilrEj and the promotion operator technique with tabulated 
vector coupling coefficients.O'L3 The vector coupling coefficients (also called Clebsch-Gordan coefficients) 
(ril^afrlril^rc) defined in 

\TxT2Tc) = l^r^) (r^abir^rc) (4) 

a . b 

and tabulated in Table A20 of Ref. (See also Table f of Ref. |77]) essentially determine how the product 
basis from two basis \Tia) and ^26) are further reduced according to the irreducible representation \Tc) 
of the point group. In this work, the latter method is adopted with some improvements to attain unique 
symmetrization. 

First of all, the overtone of the E (1/2) and F% (^3 and ^4 ) normal modes are symmetrized separately. 
For the symmetrization of the overtone of the E normal modes, the promotion operators technique with 
tabulated vector coupling coefficients has been used in Ref. |75|. The method is summarized as follows. The 
promotion operators with (E ai E b ) symmetry are the familiar creation operators which act on the basis by 

a a (E a )\n ai n b ) = y/n a + l\n a + l,n b ) 

a\{E b )\n a ,n b ) = y/nj, + l\n a ,n b + 1). (5) 

Taking the promotion operators to belong to (E a , E b ) symmetries, the symmetry of the resultant promoted 
basis function can be identified according to the vector coupling coefficients. Starting from the v=l basis 
function, one can derive all the symmetrized basis functions until v reaches a given number similar to 
climbing a v ladder. A computer code is desirable if one wants to symmetrize a high overtone. One problem 
is that spurious states with the same symmetry may arise if the promotion operator acts on every basis, 
and they are not orthogonal. For example, with promotion operators acting on the full 3^2 symmetrized 
basis (Ai + A2 + E), three nonorthogonal E states are obtained while only two are allowed for 4^2- These 
spurious states are removed by using the Gram-Schmidt orthogonalization technique. The final orthogonal 
symmetrized basis is not unique, but all possible solutions lead to symmetrized Hamiltonian matrix blocks. 

An alternative method, which yields unique symmetrization coefficients, is obtainedby finding the exact 
decomposition of the angular momentum basis function into the Cartesian basis wavefunction.Once this is 
done, the symmetrization is trivial as the angular momentum basis function is readily symmetrized according 
to the value of I modulo 3. This task seems quite complicated since the angular momentum basis functions are 
basically the associated Laeuerre polynomials and the Cartesian basis functions are essentially products of 
two Hermite polynorpials.HIj In fact, this complication can be circumvented by using the circular promotion 
operators defined byEJ 

a d = -^(4+«4) 

a l = ~^( Q a ~ ia b) ( 6 ) 

They act on the angular momentum basis by 
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and on the Cartesian basis by 



a\\n a ,n b ) = \j " 1 \n a + l,n b ) + — ' ~ \n a ,n b + 1) 



n a + 1 , 



'rift + 


1 


2 


'n 6 + 


1 



<4|n a ,n fc ) = y — K + i'"*) _ *Y — g — l nQ ' n <> + 1 ) ( 8 ) 

With the circular promotion operators acting on the v = ground state 1 0°) = |0,0) in the two basis sets, 
respectively, one obtains the relationship for the v — 1 basis 

|l ±1 > = ^(|10)±i|01)) (9) 

Repeating this operation up to a given quantum number v leads to relationship of the two basis up to v. 
The next step is to transform the angular momentum basis from the (E + , E~) to (E a , E b ) symmetries. It is 
well known that |i> ±( ) (I = 3p + 1, with p being any integer) form a basis for E^ irreducible representations. 
Following a convention which is learned from Eq. (j^), the (E a , E b ) symmetry basis functions are defined by 

E a ) = ^(\v+ l ) + \v- 1 )) = V2Re(\v+ 1 )) 
\ V W;E b ) = ^(\v +l ) - \v- 1 )) = V2Im(\v+ 1 )) (10) 
where I — 3p + 1. The A1/A2 symmetry basis functions are defined similarly by 

\vW;Ai) = ^(|« + ''l> + \v-W)) = V2Re(\v+W)) 

\vW;A 2 ) = - \v-W)) = V2Im(\v + W)) (11) 

where / = 3p. There is only A\ symmetry if I is zero. Pak et have also discussed a similar scheme to 
transform the angular momentum basis functions to the real cosine and sine form, which are members of 
the irreducible representation of the D^h group. It should be noted that the above results can be applied to 
other analogous doubly degenerated modes such as the E mode in the C$ v group and the E g mode in the Oh 
group. The above procedure can be easily implemented by-a computer code. Results up to v = 8 are listed 
in Table IV, which are obtained from a code in the MapleE3 computer algebra language. These results may 



be useful for other applications. A fortran code is also written for the same purpose and is to be combined 
with the symmctrization of the F2 normal mode to form a general symmetrization package. 

Now we turn to the symmetrization of the overtone of the F2 normal mode. This can also be achieved 
using the promotion operators together with the vector coupling coefficients. However there exist general 
rules on how a given overtone is symmetrized without resorting to promotion of lower states. These general 
rules can be derived using the projection operator methods. Roughly speaking, one obtains AE symmetry 
when the three Cartesian quantum numbers are all even or odd, and F symmetry when they are mixed even 
and odd. Detailed symmetrization coefficients have been given in Table 2 of Ref. ^0| for the overtone of the 
triply degenerated modes (Fi and F 2 modes of Tj, group) of Oh and Td group. 

Finally given the fact that the E overtone and F2 overtone have been symmetrized separately, the com- 
bination bands in the form (ni 712^714) can be symmetrized using the vector coupling coefficients, as has 
been demonstrated in Ref. |75]. In this process, no spurious states arise as the basis functions joining the 
product are from different normal modes. All the above steps have been implemented as a Fortran code 
which writes to a file the symmctrization coefficients for all the combination bands in a polyad with a given 
polyad number. This is highly desirable since one may encounter completely different polyad numbers such 
as those for CH4, CF4 and SiF4. These symmetrization coefficients are further used in constructing the 
polyad block matrix of the transformed Hamiltonian. !_. j—. 

It is interesting to compare with the basis symmetrization for the Td group in the algebraic model.E3~Ej A 
key difference is that the basis functions symmetrized there are products of ten functions based on the ten 
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internal coordinates. Due to the redundancy betsaeea the six bond-angle coordinates considerable effort was 
required to remove the so-called spurious statesSEj Their symmetrization proceduije-|is basically the same 
as the promotion operator technique with tabulated vector coupling coefficientsJZ-lO Lemus and FrankH 
have also tabul ated , and derived part of the vector coupling coefficients in the group, unaware of an earlier 
full tabulation.Ollj 



B. Extraction of spectroscopic constants 



The spectroscopic constants we want to extract from our second order transformed Hamiltonian are 
specifically the quartic anharmonicity constants Xij, Gij, the tensor constants Ty and S!y and the cubic 
Ki jk and quartic resonance constants Kijki and Kij^i- The quartic anharmonicity constants and Gij 
have their usual definitions. The tensor constants were defined by Hechttfl as coefficients of spherical tensor 
operators which cause splittings of the same symmetry components of the overtone and combination, bands of 
spherical top molecules. This procedure has been used in earlier work on the acetylene molecule£j however 
the procedure used was not described. More importantly, the procedure has not been applied to molecules 
of Td symmetry which present a special challenge. For these reasons we describe the procedure in full detail 
here. 

The apparent difficulty of the extraction is that our CVPT results for j-tha. transformed Hamiltonian are 
expressed as expansions of creation-annihilation operators in normal formOL 2 ] — the creation operators are 
put before the annihilation operators. These expansion terms need to be rearranged to obtain the traditional 
spectroscopic Hamiltonian. A simple example based on second order perturbation theory serves to illustrate 
this point. At second order the normal form expansion includes contributions, such as (aj) 2 (ai) 2 = ni(nx — 1), 
a\a\ — hi plus a constant. It is found that by rearranging the normal form terms as powers of (rii + i), 
all the linear contributions to second order of perturbation theory cancel, leaving the Hamiltonian in the 
standard spectroscopic form where only quadratic terms such as (rii + \) 2 contribute. With either choice of 
expansion the leading term will have the same coefficient which is X\\. 

A similar problem has been studied by Hodgkinson et alE3 who have suggested the same idea of extracting 
spectroscopic constants from the second order transformed Hamiltonian, although the Hamiltonian discussed 
there is in rectilinear normal coordinates and the creation-annihilation operators are not rearranged to normal 
form. Furthermore, the extraction of constants related with doubly and triply degenerated modes requires one 
more step than the above simple example of An. For example, the tensor constants T 23 and anharmonicity 
constants X 23 are determined simultaneously by 

Coef(n 2Q n 2z ) = A 23 + 8T 23 

Coef(n 2fc n 2z ) = A 23 - 8T 23 (12) 

where Coef(6) denotes coefficient of operator 0, and n 2a n 2z and n 2 b?T. 2z are the leading terms in the normal 
form Hamiltonian. These formulae are readily obtained by rewriting the 23 operator which was defined by 
Hecht [Eq. (8) of Ref. [54} as an operator whose coefficients are T 23 . Such formulae for all the anharmonicity 
and tensor constants have been given in Table 5 of Ref. [58| 

Here one has to take care of a convention problem regarding the definition of normal/symmetry coordinates 
for the V2 (E) mode of a tctrahcdral molecule. We notice that Jahn defined the symmetry coordinate 
(O20, O2&) on page 472 of Ref. ^7| (See Fig. 1 of Ref. |8| for the numbering of the four peripheral atoms). 
His definition and the vector coupling coefficients for the E mode coupling other modes (Table II of Ref. B9) 
were adopted by Hecht in Refj^ —Unfortunately the symmetry-coordinates (S^m S 1 ^) for the v% mode we 
and many other researchersE9EK'Ell adopted are those of MillsOO After some manipulations, it is found 
that the two definitions are related such that Q 2a = — >S 2 b,Q 2 b = S^o- If one uses Mills' definition, all the 
Hecht formulae regarding the v-i mode should be modified according to the above relationship. The only 
significant change occurs with the 23 operator for which a minus sign needs to be attached to Hecht 's 
definition. Only after this correction is considered, can the Table 5 of Hodgkinson et at be reconciled with 



Hecht's definition. Eq. (12) is a pertinent example. Finally, it is noted that the vector coupling coefficients 
tabulated by Tababe and SugancO and GriffithLJ agree with Mills' definition of the symmetry coordinates 

Of V2- 

The extraction of the resonance constants is straightforward in the creation-annihilation form since they 
are usually defined in this form. One just needs to take care of constructing the correctly symmetrized 
Hamiltonian terms fo r . trip resonance. To this end, several equivalent approaches can be used. Other than 



the method by HechtEJO for constructing the correctly symmetrized potential tensor operators, we choose 
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to follow an elegant and more general approach by HalonerJlS who used the tabulated vector coupling 
coefficients to construct the correctly symmetrized Hamiltonian terms from any combination of operators 
such as coordinates, momenta and creation-annihilation operators. Methane has a variety of cubic and 
quartic resonances due to its approximate 2:1 ratio of the stretch and bend frequency. Halonen has given for 
the first time all the correctly symmetrized Darling-Dennison type quartic resonance terms in Ref. |75| using 
the vector coupling methods. The resonance constants encountered in CF4 and SiF4 are limited compared 
to methane. Up to second order resonance, they are (See also Eqs. ( p^ ) - (|2l]) and ( p3| ) - (|25|)) 

(13) 
(14) 



H/hc = Ki t 22 a\(aj a + a\ b ) + h.c. 

H/hc = ^3,44 a\ x a iy a iz + al y a4 z a 4x + a\ z a ix a iy + h.c. 



for CF 4 and 



H/hc=Ki t 4,i a[(a 4x + a iy + a Ax ) + h.< 
H/hc = ifi,222 a\a 2a {a\ a - 3a 2 .,,) + h.c. 



(15) 
(16) 



for SiF4. The quartic resonance with coefficient ^1,222 is of a new type. We shall show how it is derived 
to illustrate the vector coupling method. The goal is to construct the totally symmetric (A\ symmetry) 
quartic terms from one creation operator (a\) with A\ symmetry and two annihilation operators (a 2a ,a 2 b) 
with E symmetry. Apparently there is only one such term. First the E a /Eb symmetry quadratic terms are 
constructed from products of {a 2ai a 2 b) and themselves as follows 



En 



E h 



1 2 1 2 \ 
{-a 2a + a 2b ) 



V2 
V2a 2 



(17) 



Then the A\ symmetry cubic terms are constructed from products of (a 2a , a 2 b) and the two terms of Eq. ( |l7| ) 
as follows 



;^|(- a L + a ib) ■ a 2a + V2a 2a a 2b ■ a 2b 



"2 a 2a( a 2a 



3a 2 



2b) 



(18) 



Finally the A\ symmetry quartic term is obtained as product of a\ and terms of Eq. ( |l8|) plus the Hermite 
conjugate of the product. The coefficient of the final resonance term in Eq. (16) has been redefined for 
convenience. 

An interesting point is that the anharmonicity and tensor constants extracted using the method above agree 
exactly with standard second-order perturbation theory in rectilinear normal coordinates as implemented 
in SPECTRO, provided that the same force field is used and no resonances are considered. This seems Jm 
be a surprise, since the reported CVPT results are based on curvilinear coordinates. McCoy and SiberlE3 
have explained this as follows. When using dimcnsionless normal coordinates the perturbation parameter 
in the expansion of Eq. (1) can be taken as h 1 ^ 2 . Here terms of order n in Eq. (1) are of order h!* 2+n ^ 2 . 
This is true regardless of using rectilinear or curvilinear normal coordinates. Since both sets have identical 
zero-order Hamiltonians, the energies have to be identical order by order. If one uses perturbation theory 
to transform to block-diagonal Hamiltonian, which-is subsequently diagonalized, the above argument breaks 
down, and notable differences have been found. 0'E3 In this work, a multi-resonance Hamiltonian has been 
considered for CF4 and SiF4. Therefore the spectroscopic constants are extracted from the second order 
CVPT transformed Hamiltonian, and some of them are different from those obtained using SPECTRO with 
resonances deperturbed. More details regarding this point are given with respect to CF4 and SiF4 molecules 
in the following section. 



IV. RESULTS FOR CF 4 



The choice of the polyad quantum numbers for CF4 is based on the following considerations. The harmonic 
frequencies as calculated from the ab initio force fields (u>i = 921.57 cm -1 , uj 2 = 439.91 cm -1 , W3 = 1303.01 



637.89 cm x ) suggest that three independent resonances are possible: 
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ui w 2uj 2 (19) 
w 3 ps 2lu 4 (20) 
3cj 2 ~ w 3 (21) 

The frequency differences between the above pairs are only 42 cm -1 , 27 cm" 1 and 17 cm , respectively. 
Since these small differences will go into the denominator of the expression of S^ n \ they will possibly lead 
to divergence of the perturbation theory when the coupling strength between the two pair states are not 
negligible. In practice we found that all the three resonances ought to be considered by keeping them in 
the K v to ensure good convergence when the perturbation theory is carried up to the sixth order. It should 
be noted that the resonance in Eq. ( pl| ) is actually a fourth order resonance, because the symmetry is 
unmatched between 3z^ and 1/3 at second order. Due to the three independent resonances, only one good 
quantum number remains, which should be orthogonal to the three resonance vectors (1, —2, 0, 0), (0, 0, 1, —2) 
and (0, —3, 1,0). It is then easy to find the remaining good quantum number to be 

N = 4«i + 2«2 + 6U3 + 3i>4, (22) 

which will be used in the perturbation calculations reported below. 



In Table VI, we collected all the reliable experimental band origins of the three isotopes. There are 14 
high resolution data (with resolution better .than 0.01 cm -1 ). All the other data with typical 0.1 and 1 cm -1 
resolution are from the work of Jones et al.KB There are some other low- resolution data compiled in Ref. [II], 
which we found to be unreliable and were discarded. The calculated energy levels based on the adjusted ab 



initio force field of Table [II agree well with the observations. We then further refined all five quadratic force 
constants to the observed levels with uncertainty for each level determined by the experimental precision, 
using a fast convergent second order least squares method where the Hessian is approximated as products 



of first derivatives O The refined constants are given in Table and the comparison of experimental and 



calculated energy levels is given in Table VI . The RMSD (root mean squares deviations) based on the refined 



potential is 0.179, 0.083 and 0.083 cm 1 at second, fourth and sixth order CVPT, respectively, whereas the 



corresponding RMSD based on the ab initio potential in Table [II is 0.095, 0.190 and 0.186 cm 1 . This shows 



that the ab initio potential is a rather good as an initial potential for the fitting. The following remarks can 



be made by examining the calculated results in Table VI 



(1) The CVPT calculation converges very well. The energy level changes typically on the order of 0.1 cm -1 
by comparing the fourth and sixth order results. This contrasts withJLhe case of CH4 where the bending 
overtone (4^4) varies by about three wavenumbers at the eighth order O Therefore the fourth order CVPT 
is used for the fitting. The contribution due to the V terms is small. For example, the largest contribution 
to the fundamentals is 0.057 cm -1 for (0001;F 2 ) band of 12 CF 4 using sixth order CVPT. 

(2) For 12 CF4, the resonance UJ3 ps 2w4 is the most pronounced. Even the (0010;^) fundamental has 7 
% admixture from (0002;_F 2 ). The resonance ui « 2^2 is less pronounced and the resonance 3^2 ~ can 
be ignored. The above remarks also apply for the 13 CF4 whereas all the three resonances seem to be less 
pronounced for 14 CF4. 

The pronounced resonance between W3 pa 2o>4 means-t hat p revious treatment for the ^3 ladder states where 
only the ^3 overtones are considered are inadequate E2EIE 2 ] 

(3) The calculations for the three (0020) bands using either the ab initio or fitted force field cannot 
reproduce well the observations which were obtained through a rovibrational analysis of the isolated (0020) 
bandso We cannot fully understand these discrepancies. But it is clear that the (0020) band is in Fermi 
resonance with the (0012) and (0004) bands, although the major perturbant (0012) is a little farther away. 
For example, the three (0012;^) bands are respectively at 2536.50 cm" 1 , 2541.66 cm" 1 , and 2544.45 cm" 1 . 
We have also noted that near the experimental (0020;^) band origins (2561.91 cm" 1 ), there is a band 
calculated to be 2561.04 cm" 1 with 82%(0302;F 2 ) and 14%(0310;F 2 ). However it is very unlikely that this 
band should be as strong as the (0020;^) band. Further theoretical and experimental work on this band 
system appears to be desirable. 

Finally, the vibrational spectroscopic constants for all the isotopes of CF4 and SiF4 obtained from the 



second order transformed CVPT Hamiltonian (curvilinear formalism) are given in Table VII. With these con- 
stants one can reproduce the second order CVPT calculation by constructing the block diagonal Hamilonian 
matrices. The corresponding constants in the rectilinear formalism, where different, are given in footnotes 
to said table. The differences (for the Fermi resonance constants and for such anharmonic constants as are 
affected by the Fermi interactions) illustrate the points made in the previous section. 

Discrepancies between the presently computed anharmonicity constants and the empirically derived sets 
of HeenarO mostly reflect the limitations of the latter in terms of available data. There is no indication 
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for the very high V33 = —9.1 cm -1 as suggested by Maring et a/.E3 We also note that both Heenan sets of 
harmonic frequencies differ quite substantially from our best values, and are confident that the latter are 
more reliable. 



V. RESULTS FOR SIF 4 



The choice of the polyad quantum numbers for SiF4 is based on the following considerations. The harmonic 
frequencies as calculated from the ab initio force fields (u>i = 806.3 cm -1 , lo 2 = 265.2 cm -1 , u 3 = 1044.0 
cm -1 , UJ4 = 389.3 cm -1 ) suggest that three independent resonances are possible: 

cji w 2w 4 (23) 
w\ + 0J2 ~ ^3 (24) 
0J1 w 3^2 (25) 

The differences between the above pairs are only 27 cm" 1 , 27 cm" 1 , and 11 cm" 1 , respectively. Again all 
three resonances ought to be considered to ensure good convergence when the perturbation theory is carried 
up to sixth order. It should be noted that the resonance in Eq. (E4|) is actually a third order resonance, 
because the symmetry is unmatched between v\ + v 2 (F2) and v 3 (E) at first order. Following the approach 
of the previous section the good quantum number, to be used in the perturbative calculations, is found to 
be 

N = 6vi + 2v 2 + Sv 3 + 3w 4 (26) 

which we are going to use in the perturbation calculations reported below. 

The sixteen observed band centers from Ref. pq and the values of the second, fourth and sixth order CVPT 



calculations are given in Table VIII . The following remarks can be made by examining Table VIII . 

(1) The calculated values converges to a few hundredths of a wavenumber by comparing the fourth and 
sixth order results. The convergence is better than in the case of CF4. The contribution due to the V terms 
is small. For example, the largest contribution to the fundamentals is 0.0016 cm" 1 for (0001;^) band of 
28 SiF4 using sixth order CVPT. This is smaller than that of the 12 CF4 since SiF4 is heavier. 

(2) The agreement between calculation and experiment is rather good with a RMSD of 0.59, 0.74 and 
0.73 cm -1 at second, fourth and sixth order CVPT, respectively. It should be noted that the slightly better 
agreement at the second order is fortuitous as the second order results still do not converge. On the whole, 
the CVPT results are remarkable in spite of the fact that the CVPT calculation uses the ab initio force field 
where only four diagonal quadratic force constants were optimized to the four experimental fundamentals. 
The success of the perturbation calculation is mainly due to the small anharmonicity for this relatively 
heavy-atom m olecul e (cfr. the uji — Vi values in Tables [j] and ||). 



(3) In Table VIII our calculations are compared with those of the U (2) algebraic model. El Although the 
U(2) algebraic model yields a similar RMSD (0.79 cm -1 ), nine parameters were-employed to fit the sixteen 
vibrational data. A fit of similar quality is also reported by the same authorsEZI using a boson realization 
model. However, in both models the polyad number N = v% + v 2 + v 3 + w 4 is used and the Fermi resonance 
u>i « 2^4 is considered in neither. 



(4) It can be seen from Table VIII that the resonance u>\ f=a 2^4 is most prominent. The other two 
resonances are also considered to ensure good convergence at sixth order. To neglect these two resonances 
will considerably reduce the size of the matrix associated with a given value of the quantum number, so that 
high overtone states can be relatively easier to calculate. The lack of strong resonance in v 3 and its overtone 
ensures that the isolated treatment of the v 3 overtone is a good approximation. LJ Patterson and Pine has 
done such a calculation, and find the Cartesian basis is much inferior to the angular momentum basis due 
to the near degeneracy of the .-Cartesian basis. This near degeneracy is specific to SiF4 since the relation 
G33 = 8T33 is nearly satisfiedEj However, as far as the calculation of the vibrational energy is concerned, 
both basis sets can be used. ,— , 



The partial set of empirically derived anharmonicity constants (Table VII) of McDowell et agrees 
fairly well with our computed values, while their derived harmonic frequencies agree to within their stated 
uncertainty with our own best values. 
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VI. CONCLUSIONS 



The anharmonic vibrational spectra of the all-heavy atom spherical tops CF4 and SiF4 have been treated 
by the combination of an accurate ab initio anharmonic force field and high-order canonical Van Vleck 
perturbation theory (CVPT). 

The anharmonic part of the potential energy surface is evidently very well described at the CCSD(T)/cc- 
pVTZ level for CF 4 , and at the CCSD(T)/cc-pVTZ+l level for SiF 4 , where the '+1' notation refers to 
the addition of a high-exponent d function on second-row atoms. The harmonic frequencies of CF4 exhibit 
substantial errors at the CCSD(T)/cc-pVTZ level, which disappear upon addition of diffuse functions to the 
basis set. 

The force field is subsequently slightly refined by adjusting the equilibrium geometry and the diagonal 
quadratic force constants (in this case, five parameters) such that a standard second-order rovibrational 
perturbation theoretical treatment reproduces the experimental fundamentals and ground-state rotational 
constants of the molecules. The adjustments involved are fairly minor. 

These force fields were then used as input for CVPT calculations up to sixth order inclusive. While 
agreement with experiment is fortuitously quite good at second order in CVPT, consistent convergence to 
0.1 cm -1 in the energy level is only achieved at sixth order. However, the truncation error at second order 
is much less significant than in the case of hydride systems like CH4 and SiH4. 

RMS deviation (RMSD) between computed (sixth order CVPT) and observed energy levels with the 
adjusted ab initio potentials is 0.19 cm -1 for CF4 and 0.73 cm" 1 for SiF4, the latter to some extent reflecting 
the lesser accuracy of the experimental data. In the case of CF4, improvement with experiment could be 
somewhat further improved (to RMSD=0.08 cm -1 ) by re-adjusting both diagonal and off-diagonal quadratic 
force constants to the complete set of experimental vibrational level information. Experimentation with 
refinement of additional force field parameters yielded no further improvement. 

For CF4, three resonances, Eqs. ( pa ) (p5|), were considered, of which only one (0*3 ~ 2^4) is important. 
This shows that an isolated ^3 overtone ladder model is inadequate. For SiF4, three resonances, Eqs. (pj|) 
(f2ij), were considered, of which only one ss 2^4) is important. 

An improved approach for symmetrizing combination bands in the Cartesian basis for the group is 
proposed. We also demonstrate how anharmonic spectroscopic constants Xy , Gi 3 ■, ■, Sij can be extracted 
from the second-order CVPT transformed Hamiltonian (in curvilinear internal coordinates) for molecules. 
In the absence of resonances, the results are identical to those obtained by standard second-order perturbation 
theory in rectilinear normal coordinates. Differences occur when the constants are being deperturbed for 
Fermi resonances. Accurate sets of quartic spectroscopic constants for the isotomomers of CF4 and SiF4 are 
obtained. Agreement with previously published empirically derived sets of anharmonicity constants is fairly 
good for SiF4, but less satisfactory for CF 4 . 

In order to stimulate further research on these molecules, sixth-order CVPT energy level predictions up to 
polyad number N = 24 for ^ 12 ' 13 ' 14 }CF 4 a nd f 28 ' 29 ' 3 °}SiF4 have been made available on the World Wide Web 
at the Uniform Resource Locator (URL) |ittp : //theochem. weizmann. ac . il/web/papers/cf 4sif 4 .html . 
The force fields themselves are available in machine-readable form at the same URL. 



ACKNOWLEDGMENTS 



Prof. M. S. Child is thanked for helpful discussions on the basis set symmetrization. JM is a Yigal Allon 
Fellow and the incumbent of the Helen and Milton A. Kimmclman Career Development Chair, acknowledges 
support from the Minerva Foundation, Munich, Germany, and would like to thank Dr. Timothy J. Lee 
(NASA Ames Research Center) for helpful discussions and Prof. Jean Demaison (Universite de Lille I) for 
a preprint of Ref. S XGW thanks the Royal Society KC Wong Fellowship for support. 



1 Kirk-Othmer encyclopedia of chemical technology, 4th Ed., R. E. Kirk, D. F. Othmer, J. I. Kroschwitz, Eds. (Wiley, 
New York, 1998) 

2 H.-G. Mack and H. Oberhammer, J. Chem. Phys. 87, 2158 (1987) and references therein 



11 



3 J.-H. Kim, S.-H. Seo, S.-M. Yun, H.-Y. Chang, K.-M. Lee, and C.-K. Choi, J. Electrochem. Soc. 143, 2990 (1996) 
and references therein. 

4 J. Ding and N. Hershkowitz, Appl. Phys. Lett. 68, 1619 (1996) 

5 P. G. M. Sebel, L. J. F. Hermans, and H. C. W. Beijerinck, J. Vac. Sci. Tech. A 17, 755 (1999); see also T. Lagarde, 
J. Pelletier, and Y. Arnal, J. Vac. Sci. Tech. A 17, 118 (1999) 

6 P. Francis, C. Chaffin, A. Maciejewski, and C. Oppenheimer, Geophys. Res. Lett. 23, 249 (1996) 

7 K. Tanaka, S. Isomura, H. Kaetsu, Y. Yatsurugi, M. Hashimoto, K. Togashi, and S. Arai, Bull. Chem. Soc. Japan 
69, 493 (1996) 

8 B. Monostori and A. Weber, J. Chem. Phys. 33, 1867 (1960); A. Maki, E. K. Plyler, and R. Thibault, J. Chem. 
Phys. 37, 1899 (1960); P. J. H. Woltz and A. H. Nielsen, J. Chem. Phys. 20, 307 (1952); 

9 L. H. Jones, B. J. Krohn, and R. C. Kennedy, J. Mol. Spectrosc. 70, 288 (1978) 

10 L. H. Jones, C. Kennedy, and S. Ekberg, J. Chem. Phys. 69, 833 (1978). 

11 A. C. Jeannotte II, D. Legler, and J. Overend, Spectrochim. Acta A 29, 1915 (1973) 

12 J. J. Tiee and C. Wittig, Appl. Phys. Lett. 30, 420 (1977); J. Appl. Phys. 49, 61 (1978) 

13 P. Esherick, A. Owyoung, and C. W. Patterson, J. Molec. Spectrosc. 86, 250 (1981). 

14 J. E. Lolck, J. Raman Spectrosc. 11, 294 (1981). 

lj A. Tabyaoui, B. Lavorel, R. Saint-Loup, and M. Rotger, J. Raman Spectrosc. 25, 255 (1994). 

16 C. W. Patterson, R. S. McDowell, N. G. Nereson, R. F. Begley, H. W. Galbraith, and B. J. Krohn, J. Molec. 
Spectrosc. 80, 71 (1980). 

17 G. Poussigue, G. Tarrago, and A. Valentin, J. Phys. B. 47, 1155 (1986) 

18 G. Tarrago, G. Poussigue, M. Dang-Nhu, and J. Kauppinen, J. Mol. Spectrosc. 86, 232 (1981) 

19 R. S. McDowell, M. J. Reisfeld, H. W. Galbraith, B. J. Krohn, H. Flicker, R. C. Kennedy, J. P. Aldridge, and 
N. G. Nereson, J. Molec. Spectrosc. 83, 440 (1980). 

20 A. S. Pine, J. Molec. Spectrosc. 96, 395 (1982) 

21 M. Dang-Nhu, G. Graner, and G. Guelachvili, J. Mol. Spectrosc. 85, 358 (1981) 

22 M. Takami, J. Chem. Phys. 73, 2665 (1980) 

23 T. Gabard, A. Nikitin, J. P. Champion, G. Pierre, and A. S. Pine, J. Molec. Spectrosc. 170, 431 (1995). 

24 T. Gabard, L. Pierre, and M. Takami, Mol. Phys. 85, 735 (1995) 

25 J. L. Duncan and I. M. Mills, Spectrochim. Acta 20, 1089 (1964) 

26 A. A. Chalmers and D. C. McKean, Spectrochim. Acta 22, 251 (1966) 

27 A. C. Jeannotte II, C. Marcott, and J. Overend, J. Chem. Phys. 65, 2076 (1978) 

28 S. Brodersen, J. Mol. Spectrosc. 145, 331 (1991) 

29 S. G. Larsen and S. Brodersen, J. Molec. Spectrosc. 157, 220 (1993). 

30 V. Boujut, F. Michelot, and C. Leroy, Mol. Phys. 93, 879 (1998). 

31 V. G. Tyuterev, Y. L. Babikov, S. A. Tashkun, V. I. Perevalov, A. Nikitin, J. P. Champion^. Wenger, G. Pierre, 
J. C. Hilico, and M. Loete, J. Quant. Spec. Radiat. Transfer 52, 459 (1994) based on Ref. |2J 

32 W. Maring, J. P. Toennies, R. G. Wang, and H. B. Levene, J. Chem. Phys. 103, 1333 (1995). 

33 R. K. Heenan, Ph.D. Thesis (U. of Reading, UK, 1979) 

34 K. T. Hecht, J. Molec. Spectrosc. 5, 355 (1960). 

35 R. S. McDowell, M. J. Reisfeld, C. W. Patterson, B. J. Krohn, M. C. Vasquez, and G. A. Laguna, J. Chem. Phys. 
77, 4337 (1982). 

36 X.-W. Hou, S.-H. Dong, M. Xie, and Z.-Q. Ma, Chem. Phys. Lett. 283, 174 (1998). 
37 X.-W. Hou, M. Xie, S.-H. Dong, and Z.-Q. Ma, Ann. Phys. (NY) 203, 340 (1998). 

38 F. Iachello, Chem. Phys. Lett. 78, 581 (1981); F. Iachello and R. D. Levine, J. Chem. Phys. 77, 3046 (1982) 

39 C. W. Patterson and A. S. Pine, J. Mol. Spectrosc. 96, 404 (1982) 

40 M. Takami and H. Kuze, J. Chem. Phys. 78, 2204 (1983). 

41 L. Jorissen, H. Prinz, W. A. Kreiner, C. Wenger, G. Pierre, G. Magerl, W. Schupita, Can. J. Phys. 67, 532 (1989). 

42 J. Breidung, J. Demaison, L. Margules, and W. Thiel, Chem. Phys. Lett., submitted ("Equilibrium Structure of 
SiF4"); preprint communicated to authors 

43 J. M. L. Martin and P. R. Taylor, J. Phys. Chem. A 103, 4427 (1999). 

44 G. K. Johnson, J. Chem. Thermodyn. 18, 801 (1986) 

45 T. J. Lee, J. M. L. Martin, and P. R. Taylor, J. Chem. Phys. 102, 254 (1995) 

46 S. Carter, H. M. Shnider, and J. M. Bowman, J. Chem. Phys. 110, 8417 (1999) 

47 E. Venuti, L. Halonen, and R. G. Delia Valle, J. Chem. Phys. 110, 7339 (1999). 

48 X.-G. Wang and E. L. Sibert, J. Chem. Phys. Ill, 4510 (1999). 

49 H.-J. Werner and P. J. Knowles, MOLPRO 97.3, a package of ab initio programs, with contributions from J. 
Almlof, R. D. Amos, A. Berning, D. L. Cooper, M. J. O. Deegan, A. J. Dobbyn, F. Eckert, S. T. Elbert, C. 
Hampel, R. Lindh, A. W. Lloyd, W. Meyer, A. Nicklass, K. A. Peterson, R. M. Pitzer, A. J. Stone, P. R. Taylor, 
M. E. Mura, P. Pulay, M. Schiitz, H. Stoll and T. Thorsteinsson. 

50 G. D. Purvis III and R. J. Bartlett, J. Chem. Phys. 76, 1910 (1982) 



12 



51 K. Raghavachari, G. W. Trucks, J. A. Pople, and M. Head-Gordon, Chem. Phys. Lett. 157, 479 (1989) 

52 P. J. Knowles, C. Hampel, and H. J. Werner, J. Chem. Phys. 99, 5219 (1993) 

53 T. J. Lee and G. E. Scuseria, in Quantum mechanical electronic structure calculations with chemical accuracy, Ed. 
S. R. Langhoff (Kluwer Academic Publishers, Dordrecht, The Netherlands, 1995). 

54 T. J. Lee and P. R. Taylor, Int. J. Quantum Chem. Symp. 23, 199 (1989) 

55 D. Papousek and M. R. Aliev, Molecular Vibrational- Rotational Spectra (Elsevier Scientific Publishing Company,: 
New York 1982). 

56 A. Willetts, J. F. Gaw, W. H. Green Jr., and N. C. Handy, SPECTRO 1.0, a second-order rovibrational perturbation 
theory program (University Chemical Laboratory, Cambridge, UK, 1989) 

57 J. F. Gaw, A. Willetts, W. H. Green, and N. C. Handy, in Advances in molecular vibrations and collision dynamics 
(ed. J. M. Bowman), JAI Press, Greenwich, CT, 1990. 

58 D. P. Hodgkinson, R. K. Heenan, A. R. Hoy, and A. G. Robiette, Mol. Phys. 48, 193 (1983). 

59 This is fairly trivially done by taking apart the compound fractions in Eq.(3.20) of Ref. |5£|, then inserting a test that 
skips the appropriate one among the four affected terms if the respective denominator corresponds to the resonant 
interaction. Alternatively, the reciprocal denominators may be precomputed and stored in a 3-dimensional array, 
and any that correspond to resonant interactions set to zero before the actual spectroscopy routines are entered, 
similar to J. M. L. Martin and P. R. Taylor, Spectrochim. Acta A 53, 1039 (1997). 

60 T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989) 

61 D. E. Woon and T. H. Dunning Jr., J. Chem. Phys. 98, 1358 (1993). 

62 J. M. L. Martin and O. Uzan, Chem. Phys. Lett. 282, 16 (1998) 

63 J. M. L. Martin, J. Chem. Phys. 108, 2791 (1998) 

64 J. M. L. Martin, Spectrochim. Acta A 55, 709 (1999) (special issue "Theoretical spectroscopy: state of the science") 

65 R. A. Kendall, T. H. Dunning Jr., and R. J. Harrison, J. Chem. Phys. 96, 6796 (1992). 
J. M. L. Martin and P. R. Taylor, Chem. Phys. Lett. 225, 473 (1994). 
J. M. L. Martin, J. P. Francois, and R. Gijbels, J. Phys. Chem. 98, 11394 (1994). 
J. M. L. Martin, Chem. Phys. Lett. 292, 411 (1998) 

J. M. L. Martin, K. K. Baldridge, and T. J. Lee, Mol. Phys. xx, yyy (1999) 
H. M. Pickett, J. Chem. Phys. 56, 1715 (1972). 
E. L. Sibert, J. Chem. Phys. 88, 4378 (1988). 
E. L. Sibert, Comp. Phys. Comm. 51, 149 (1988). 

R. G. Delia Valle, L. Halonen, and E. Venuti, J. Comput. Chem., (in press). 
L. Halonen and M. S. Child, Comp. Phys. Comm. 51, 173 (1988). 
L. Halonen, J. Chem. Phys. 106, 831 (1997). 

J. S. Griffith, The Theory of Transition- Metal Ions, Cambridge University Press, Cambridge, (1961), Table A20 
of Appendix 2. 

Y. Tanabe and S. Sugano, J. Phys. Soc. Jap. 9, 753 (1954). 
78 L. Pauling and E. B.Wilson, Introduction to Quantum Mechanics with Applications to Chemistry (Dover Publica- 
tions, Inc.: New York 1935). 
W. H. Shaffer, Rev. Mod. Phys. 16, 245 (1944). 

C. Cohen- Tannoudji, B. Diu, and F. Laloe, Quantum Mechanics, vol 1, pp 727-741, Wiley, New York, (1977). 
Y. Pak, E. L. Sibert III, and R. C. Woods, J. Chem. Phys. 107, 1717 (1997). 
MAPLE V Release 4, Waterloo Maple Software, Waterloo, Ontario, 1981-1994. 
R. Lemus and A. Frank, J. Chem. Phys. 101, 8321 (1994). 
Z.-Q. Ma, X.-W. Hou, and M. Xie, Phys. Rev. A 53, 2173 (1996). 
85 A. Frank, R. Lemus, F. Perez-Bernal, and R. Bijker, J. Mol. Spectrosc. 196, 329 (1999). 

A. B. McCoy and E. L. Sibert III, J. Chem. Phys. 105, 459 (1996); E. L. Sibert III and A. B. McCoy, J. Chem. 
Phys. 105, 469 (1996). 

H. A. Jahn, Proc. Roy. Soc. A168, 469 (1938). 
H. A. Jahn, Ann. Phys. Lpz. 23, 529 (1935). 

H. A. Jahn, Proc. Roy. Soc. A168, 495 (1938). 

90 D. L. Gray and A. G. Robiette, Mol. Phys. 37, 1901 (1979). 

I. M. Mills, Mol. Phys. 1, 107 (1958), Fig. 1 of this work is inconsistent with its definition of symmetry coordinates. 
A. G. Robiette, D. L. Gray, and F. W. Birss, Mol. Phys. 32, 1591 (1976). 
A. B. McCoy and E. L. Sibert, Mol. Phys. 77, 697 (1992). 
A. B. McCoy and E. L. Sibert, J. Chem. Phys. 95, 3476 (1991). 

95 A. B. McCoy, D. C. Burleigh, and E. L. Sibert, J. Chem. Phys. 95, 7449 (1991). 
K. L. Mardis and E. L. Sibert, J. Molec. Spectrosc. 187, 167 (1998). 



77 



!)(> 



13 



TABLE I. Computed (CCSD(T)) and observed bond distances (A), harmonic and fundamental frequencies (cm l ) 



for 12 CF 4 





cc 


pVTZ unadj 




aug-cc-pVTZ 


cc- 


pVTZ adjusted 


Expt.S 






r e 




r x 


r e 


r e 


r 


r z 


r e r 






1.31919 


1.32310 


1.32389 


1.32112 


1.31526 


1.31925 


1.32004 


1.3151(17) a 1.319247(1)" 


i 


Uli 


Vi 


Uli — Vi 


Uli 


Uli 


Vi 


Uli — Vi 


Uli Vi 






1 


922.80 


910.75 


12.06 


915.2 


921.57 


909.07 


12.50 


909.0720(1) 






2 


440.00 


435.59 


4.40 


435.2 


439.91 


435.40 


4.51 


435.399(10) 


14 




3 


1322.25 


1303.15 


19.10 


1301.3 


1303.01 


1283.66 


19.35 


1283.66429(12.) 


-A 


4 


638.81 


632.22 


6.59 


630.4 


637.89 


631.06 


6.84 


631.05890(13 


\|24| 



Rcf. § 



b From B o =0. 19118709(32) cm -1 !! 

TABLE II. Computed (CCSD(T)) and observed bond distances (A), harmonic and fundamental frequencies (cm -1 ) 



for 28 SiF 4 





cc- 


aug-cc- 




















pVTZ 


pVTZ+1 


ce- 


pVTZ+1 unadj. 


cc-pVTZ+1 adjusted 


Expt. 






r e 


r E 


re 


ro 


r z 


r e 


ro 


r z 


r e 


ro 




1.56949 


1.56332 


1.56134 


1.56368 


1.56453 


1.55182 


1.55404 


1.55489 


1.5524(8)" 


1.55404" 


1 


U>i 


U>i 


Uli 


Vi 


Uli - Vi 


Uli 


Vi 


Uli - Vi 


U>i 


Vi 


1 


794.9 


794.1 


797.86 


792.19 


5.67 


806.10 


800.60 


5.50 


807.1(12) 


800.6 


2 


259.8 


258.8 


263.18 


262.13 


1.05 


265.20 


264.20 


1.00 


267(3) 


264.2 


3 


1036.5 


1029.8 


1037.49 


1024.31 


13.18 


1044.04 


1031.40 


12.64 


1044.2(12) 


1031.3968 


4 


384.3 


382.7 


387.61 


386.67 


0.94 


389.31 


388.44 


0.87 


389.8(9) 


388.4448 



Rcf. y. 



6 From Ref. || B o =0. 137780439(92) cm -1 . In older workHS, £ =0.13676(3) cm" 1 and hence 
ro = l. 55982(17) A. 



TABLE III. Quadratic, cubic and quartic force constants (aJ/A m radian") for S1F4 and CF4 





SiF 4 


CF 4 




SiF 4 


CF 4 




SiF 4 


CF 4 


F n 


7.27355 


9.50679 


F22 


0.63194 


1.24913 


F44 


6.40971 


6.13519 


F 74 


-0.34328 


-1.03126 


F77 


1.06468 


1.83767 


Fm 


-20.29450 


-32.46123 


F211 


-0.69851 


-2.19918 


F441 


-18.73201 


-23.71962 


F741 


0.56578 


2.57073 


F771 


-0.97378 


-3.04180 


-P222 


-0.53376 


-1.11062 


i<662 


-0.60404 


-2.24497 


Fg62 


0.61671 


2.05526 


-F992 


-1.01287 


-2.25675 


F 6 54 


-18.17256 


-18.94290 


F954 


-0.17620 


-0.38638 


-F984 


-0.01658 


0.04702 


Fgsr 


0.86987 


1.30646 




50.47669 


89.63168 


-F2211 


1.26945 


4.44699 


F4411 


48.94593 


71.68647 


-F7411 


-0.80741 


-5.46450 


-F7711 


1.53293 


5.89090 


F2221 


0.76478 


2.29325 


-^6621 


1.16338 


5.73913 


F9621 


-1.08209 


-4.08511 


J9921 


1.46146 


4.33621 


■^6541 


48.23448 


60.63713 


-F9541 


0.40465 


1.24159 


^9841 


0.15861 


0.26676 


Fga71 


-0.63502 


-1.85550 


F2222 


0.76756 


2.44329 


^6622 


0.76815 


3.16961 


^6633 


-0.46055 


-1.86075 


Fg@22 


-1.12641 


-3.26978 


^9633 


0.03132 


0.14651 


-F9922 


2.11522 


5.37994 


F9933 


0.39929 


0.77598 


F9542 


0.43597 


1.12997 


-^8762 


0.01923 


0.04530 


-F4444 


50.01306 


72.17993 


-F5544 


48.79060 


59.28204 


F7444 


-0.85975 


-6.35131 


-F8544 


0.39771 


4.59345 


-FV744 


1.37066 


5.70508 


Fb754 


0.02280 


-0.04550 


-F&844 


-0.53371 


-1.07258 


F7774 


-2.21128 


-6.91109 


F8874 


-0.14830 


-0.36797 


-F7777 


4.58572 


11.28388 


i*8877 


2.18342 


3.67611 
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TABLE IV. The relationships of the angular momentum basis \v l ) and the Cartesian basis \n x n y ) for the overtone 
of the doubly degenerated modes (1^2). See Eqs. ( p"o| ) and ( pj| ) for definition of A1/A2 and E a /Eb symmetry 
wavefunction. 
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E b) 



,72 
~ 2 
_V5 
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+ 1 

~ 2 

_vj 
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I V3 
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TABLE V. Comparison of ab initio a and fitted quadratic force constants for CF 4 . (Units are consistent with aJ, 
A, and radian.). 



Force Constant 


ab initio 


ab initio 


Fitted" 


a d 




raw" 


adjusted 6 






r e 


1.31919 


1.31526 


1.31526 




Fn 


9.54112 


9.50679 


9.50711 


0.0% 


F22 


1.25733 


1.24913 


1.12477 


-10.0% 


F33 


6.32376 


6.13519 


6.18043 


+0.7% 


F 34 


-1.03160 


-1.03126 


-1.04559 


-1.4% 


F44 


1.84886 


1.83767 


1.82662 


-0.6% 



a Unadjusted CCSD(T)/cc-pVTZ values 



b Bond distance and diagonal quadratics adjusted to reproduce experimental tq and vi 

c bond distance held constants; all quadratic force constants refined in fit againt complete experimental data 

set. 

d Relative deviation of the fitted force constants from the adjusted ab initio force constants. 
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TABLE VI. Comparison of experimental and CVPT a band origins (cm ) for CF4. The fit quadratic force 
constants of Table M plus the ab initio cubic and quartic force contants of Table pTll are used for the CVPT calculation. 



N a 


Sym 


Obs b 


Uncertainty Rcf. Obs-B(2) 


Obs-£(4) 


Obs-£(6) 


v 1 v 2 v 3 v 4 (cl) d 


"CF 4 






2 


E 


435.399 


1 


i~ 


.227 


.100 


.103 


0100(100%) 




3 


F 2 


631.0593 


1 


771 
Ml 


.127 


-.039 


-.031 


0001(100%) 




4 


Al 


867.90588 


1 


TV 


.476 


-.012 


-.003 


0200(97%) 


1000(2%) 


4 


Ai 


909.0720 


1 


Li 


-.151 


-.044 


-.045 


1000(97%) 


0200(2%) 


5 


Fi 


1066.6977 


1 


TT 


.300 


.052 


.065 


0101(100%) 




5 


F 2 


1066.1220 


1 


TT 
ill 


.320 


-.135 


-.112 


0101(100%) 




6 


A\ 


1261.809 


10 


Ml 


.303 


-.184 


-.148 


0002(90% ) 




6 


E 


1262.112 


1 


21 


.219 


-.244 


-.221 


0002(100%) 




6 


F 2 


1260.430 


1 


31 


.186 


.173 


.188 


0002(93%) 


0010(7%) 


6 


F 2 


1283.720 


1 


23 


.361 


-.005 


-.002 


0010(93%) 


0002(7%) 


7 


F 2 


1539.3 


10 


n 


-.198 


-.251 


-.241 


1001(97%) 


0201(3%) 


8 


F-2 


1715.8 


10 


n 


-.225 


-.963 


-.960 


0110(91%) 


0102(9%) 


9 


F 2 


1889. 6" 




m 


.623 


.887 


.912 


0003(84%) 


0011(16%) 


9 


F 2 












1893.42 


0003(100%) 




9 


F 2 


1913.2 


10 




2.031 


1.050 


1.100 


0011(84%) 


0003(16%) 


10 


F 2 


2168.5 


10 


n 


.674 


.604 


.621 


1002(89%) 


1010(9%) 


10 


F 2 


2186.1 


10 


tt 


-.602 


-.636 


-.632 


1010(86%) 


1002(9%) 


11 


F 2 


2445.59644 


1 




-.340 


.025 


.051 


2001(94%) 


1201(5%) 


12 


A! 


2553.24(858) 




2- 


6.225 


3.072 


3.215 


0020(50%) 


0012(41%) 


12 


E 


2570.013 






1.991 


1.652 


1.654 


0020(93%) 


0012(7%) 


12 


F 2 


2561.9124 




2: 


-2.127 


-3.061 


-3.067 


0020(86%) 


0012(13%) 


1J CF 4 






3 


F 2 


629.2868 


1 


Lfl 


.270 


.110 


.118 


0001(100%) 




5 


F 2 


1064.39 


10 


n 


.408 


-.033 


-.011 


0101(100%) 




6 


F 2 


1241.7 


10 


u 


-.251 


.014 


.009 


0010(91%) 


0002(9%) 


6 


F 2 


1259.75 


10 


n 


.775 


.138 


.160 


0002(91%) 


0010(9%) 


7 


F 2 


1537.4 


10 


n 


-.078 


-.116 


-.107 


1001(97%) 


0201(3%) 


8 


F 2 


1674.7 


10 


u 


-.099 


-.081 


-.087 


0110(93%) 


0102(7%) 


9 


F 2 


1867. 


100 


n 


1.074 


1.631 


1.630 


0011(89%) 


0003(11%) 


9 


F 2 


1888. " 




TT 


1.262 


.279 


.353 


0003(100%) 




9 


F 2 












1889.37 


0003(89%) 


0011(11%) 


10 


F 2 


2145. 


100 


u. 


-.698 


-.518 


-.530 


1010(92%) 


1002(5%) 


10 


F 2 


2166.3 


10 


m 


.017 


-.232 


-.200 


1002(92%) 


1010(5%) 


11 


F 2 


2443.3 


10 


u 


-.512 


.011 


.048 


2001(94%) 


1201(6%) 


12 


F 2 


2477.5 


10 


IT 


-3.086 


-1.556 


-1.565 


0020(71%) 


0012(26%) 


14 CF 4 






3 


F 2 


627.3490 


1 


lfl 


.275 


.121 


.128 


0001(100%) 




5 


F 2 


1062.57 


10 


m 


.449 


.019 


.041 


0101(100%) 




6 


F 2 


1208.7 


10 


n 


-.107 


-.055 


-.059 


0010(99%) 




6 


F 2 


1254.95 


10 


n 


.800 


.387 


.407 


0002(99%) 




7 


F 2 


1535.3 


10 


m 


-.204 


-.230 


-.221 


1001(97%) 


0201(3%) 


8 


F 2 


1641.6 


10 


nj 


-.038 


-.239 


-.246 


0110(99%) 




9 


F 2 


1833.4 


10 


n 


1.846 


2.023 


2.030 


0011(98%) 


0003(2%) 


9 


F 2 


1881. 4 e 




m 


.717 


-.212 


-.136 


0003(99%) 


0011(1%) 


9 


F 2 












1882.11 


0003(99%) 




10 


F 2 


2073. 7 f 


10 


if 


1.721 


1.109 


1.098 


0210(97%) 


1010(3%) 


10 


F 2 












2075.77 


0210(99%) 




10 


F 2 


2112. 


100 


M 


-.683 


-.538 


-.547 


1010(97%) 


0210(3%) 


10 


F 2 


2161.9 


10 


n 


.398 


.218 


.246 


1002(96%) 


0202(3%) 


12 


F 2 


2412. 


100 


if 


-3.208 


-2.645 


-2.638 


0020(97%) 


0012(2%) 


RMSD 










0.0179 


0.083 


0.083 






a N = 


4vi + 


2V2 + 6«3 + 3V4 


This polyad number results from three independent 


resonances: 





3lu2 ~ W3, and LU3 « 2w4. 

The last figure is significant unless uncertainty in parenthesis is given otherwise. 
c Uncertainty used in the fit corresponds approximately to the experimental precison. The states excluded 
from the fit are 3i/4 and 2i/ 2 + v% whose multiple F 2 components are unresolved, and the 2f 3 bands of 12 CF 4 , 



which differs from experimental results. 

d The largest two components in terms of V1V2V3V4 and their percentage (c|) in the basis, based on the sixth 
order calculation. Only the components with the percentage larger than 1% are listed. 
e There are two (0003) F 2 states. 
f There are two (0210) F 2 states. 
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TABLE VII. Computed (second order transformed CVPT Hamiltonian) and experimentally derived vibrational 
spectroscopic constants of S1F4 and CF4 (cm -1 ). 



Constants 



^CF 4 



II 



3 CF 4 



4 CF 4 



B SiF 4 



28 Si 
Ref. 



u SiF 4 



u SiF 4 



U3\ 


921.596 


910.52 


913.25 


921.596 


921.596 


806.101 


807.1(12j 


806.101 


806.101 


OJ2 


439.665 


437.15 


435.37 


439.665 


439.665 


265.199 


267(3) 


265.199 


265.199 


OJ3 


1302.510 


1292.96 


1302.26 


1262.587 


1227.629 


1044.044 


1044.2(12) 


1034.820 


1026.177 




637.681 


634.35 


635.83 


635.548 


633.406 


389.306 


389.8(9) 


387.793 


386.320 


X11 


-1.109 


-0.279 


-0.275 


-1.109 


-1.109 


-.669 


-0.57(5) 


-.669 


-.669 




-.911" 


0.589 


0.474 


-.911 


-.911 


-.317 


-0.6(11) 


-.317 


-.317 


X13 


-6.253 


-4.896 


-4.954 


-5.908 


-5.668 


-3.573 


-3.8(3) 


-3.503 


-3.444 


X14 


-.675 


0.608 


0.712 


-.928 


-1.087 


.109° 


+0.64(7) 


.056 


.013 


X22 


-.339° 


-0.284 


-0.254 


-.339 


-.339 


-.168 


0.0(5) 


-.168 


-.168 


X23 


-2.372 


-0.544 


-0.543 


-2.315 


-2.264 


-.989 


-1.5(10) 


-.979 


-.968 


X24 


-.004 


-0.283 


-0.273 


.003 


.007 


.506 


0.0(5) 


.503 


.500 


X33 


-4.157 


-5.642 


-5.536 


-3.917 


-3.713 


-2.845 


-3.0058(7) 


-2.779 


-2.717 


X34 


-3.960° 


-1.407 


-1.423 


-3.756 


-3.594 


-1.386 


-0.5(4) 


-1.350 


-1.318 


X44 


-.129° 


-0.606 


-0.604 


-.122 


-.114 


.165° 


-0.22(10) 


.163 


.161 


G22 


.542° 


0.286 


0.257 


.542 


.542 


.391 


[0] 


.391 


.391 


G33 


3.833 


4.317 


4.241 


3.618 


3.432 


1.799 


1.7828(6) 


1.746 


1.698 


G34 


-.518 


-1.270 


-1.320 


-.442 


-.401 


-.414 


f» -0.5 


-.377 


-.345 


G44 


.025° 


0.088 


0.084 


.025 


.025 


.000° 


-0.05(15) 


.000 


.001 


T23 


-.039 


0.045 


0.044 


-.037 


-.035 


-.101 


[0] 


-.098 


-.096 


T24 


-.037 


-0.006 


-0.005 


-.038 


-.039 


.005 


[0] 


.004 


.003 


T33 


.221 


0.683 


0.670 


.221 


.222 


.196 


0.20292(10) 


.195 


.195 


T34 


.080° 


0.391 


0.394 


.065 


.051 


-.133 


R3 -0.5 


-.131 


-.129 


T44 


.020° 


0.091 


0.090 


.021 


.021 


.014 


[0] 


.014 


.013 


S34 


.157° 






.211 


.243 


.019 


[0] 


.027 


.033 


if 1,44 















2.493° 




2.465 


2.438 


Kl,22 


-3.418° 






-3.418 


-3.418 













^3,44 


-5.223° 






-4.524 


-3.941 













^1,222 















0.063 




0.063 


0.063 



Spectroscopic constants are identical between curvilinear and rectilinear formalisms except for those affected by resonances: 

a Due to oji ps 2uj2 resonance, in rectilinear formalism: Xx2 — —0.606; X22 — —0.416; G22 — 0.619; and ^1,22 — —2.908 cm -1 . 

b Due to LJ3 PS 20J4 resonance, in rectilinear formalism: X34 — —3.916; X44 — —0.142; G34 — —0.550; G44 — 0.018; T34 — 0.067; 

T 44 = 0.024; S 34 = 0.161; and K 3 ,44 = 5.051 cm" 1 . 

c Due to resonance u>i S3 2lj4, in rectilinear formalism: X14 — —0.231; X44 — 0.250; G44 — —0.085; and K± t 44 — 2.699 cm" 1 . 
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TABLE VIII. Comparison of experimental, CVPT a and U(2) algebraic model band origins (cm 1 ) for SiF4. The 
ab initio force contants of Table pTA are used for the CVPT calculation. 



N" 


Sym 


Obs 6 


Obs-F(2) 


Obs-£(4) 


Obs-F(6) 


Obs-[/(2) c 


v t u 2 u i v A (cf) d 




2 


E 


264.2(10) 


.00 


-.03 


-.03 


-.21 


0100(100%) 




:s 


F 2 


388.4448(2) 


-.01 


-.03 


-.03 


-.41 


0001(100%) 




6 


Ai 


800.6(3) 


-.24 


-.30 


-.30 


1.03 


1000(93%) 


0002(6%) 


6 


F 2 


776.3(5) 


-.81 


-.84 


-.83 


.97 


0002(100%) 




8 


F 2 


1031.3968(3) 


-.64 


-.76 


-.75 


1.72 


0010(100%) 




8 


E 


1064.2(4) 


.00 


-.06 


-.06 


.01 


1100(92%) 


0102(7%) 


9 


F 2 


1164.2(2) 


.48 


.48 


.51 


.03 


0003(90%) 


1001(10%) 


9 


F 2 


1189.7(3) 


-.70 


-.83 


-.82 


1.07 


1001(90%) 


0003(10%) 


10 


E, 


1294.05(10) 


.25 


.11 


.11 


.15 


0110(100%) 




11 


E 2 


1418.75(10) 


.16 


.08 


.08 


.22 


0011(100%) 




14 


E, 


1804.5(1)" 


1.05 


.70 


.72 


-.21 


0012(97%) 


1010(2%) 


14 


F 2 




1805.34 


1805.49 


1805.48 




0012(96%) 


1010(4%) 


14 


Ei 




1807.20 


1807.38 


1807.37 




0012(100%) 




14 


Ei 


1828.17(2) 


-.56 


-.65 


-.65 


-.57 


1010(93%) 


0012(7%) 


16 


Ei 


2059.1(3) 


-.04 


-.17 


-.18 


1.09 


0020(100%) 




20 


Ei 


2602.55(10) / 


-.14 


-.71 


-.67 


-1.24 


1012(70%) 


2010(12%) 


20 


Ei 




2604.97 


2605.28 


2605.09 




1012(81%) 


0014(9%) 


20 


Ei 




2606.96 


2607.20 


2607.04 




1012(86%) 


0014(14%) 


20 


Ei 


2623.8(1) 


-1.30 


-1.47 


-1.43 


.12 


2010(83%) 


1012(15%) 


24 


F 2 


3068.5(1) 


-.74 


-1.71 


-1.68 


-.63 


0030(100%) 




RMSD 






0.59 


0.74 


0.73 


0.79 







a N — 6v± + 2v2 + 8V3 + 3i>4. This polyad number results from three independent resonances: u)\ + LU2 ~ UI3, 
u>i w 2^4, and w% ~ 3w2- 

b The observations are from Ref. ^5|. Standard deviations given in parentheses. 
c The U(2) algebraic model calculations are from Ref. [36|. 

d The largest two components in terms of i>ii>2^3i'4 and their percentage (c?) in the basis, based on the sixth 
order calculation. Only the components with the percentage larger than 1% are listed. 
e There are three (0012) F 2 states. 
/ There are three (1012) F 2 states. 
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